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The small-scale dynamo provides a highly efficient mechanism for the conversion of turbulent into 
magnetic energy. In astrophysical environments, such turbulence often occurs at high Mach numbers, 
implying steep slopes in the turbulent spectra. It is thus a central question whether the small-scale 
dynamo can amplify magnetic fields in the interstellar or intergalactic media, where such Mach 
numbers occur. To address this long-standing issue, we employ the Kazantsev model for turbulent 
magnetic field amplification, systematically exploring the effect of different turbulent slopes, as 
expected for Kolmogorov, Burgers, the Larson laws and results derived from numerical simulations. 
With the framework employed here, we give the first solution encompassing the complete range of 
magnetic Prandtl numbers, including Pm <C 1, Pm ~ 1 and Pm S> 1. 

We derive scaling laws of the growth rate as a function of hydrodynamic and magnetic Reynolds 
number for Pm <^ 1 and Pm S> 1 for all types of turbulence. A central result concerns the 
regime of Pm ~ 1, where the magnetic field amplification rate increases rapidly as a function of 
Pm. This phenomenon occurs for all types of turbulence we explored. We further find that the 
dynamo growth rate can be decreased by a few orders of magnitude for turbulence spectra steeper 
than Kolmogorov. We calculate the critical magnetic Reynolds number Rmc for magnetic field 
amplification, which is highest for the Burgers case. As expected, our calculation shows a linear 
behaviour of the amplification rate close to the threshold proportional to (Rm — Rmc). Based on the 
Kazantsev model, we therefore expect the existence of the small-scale dynamo for any given value 
of Pm, as long as the magnetic Reynolds number is above the critical threshold. 



I. INTRODUCTION 



The conversion of kinetic energy into magnetic energy, 
the so-called dynamo action, plays an important role in 
a wide range of astrophysical applications. These include 
the amplification of magnetic fields in the surface of the 
Sun [ll,[2, the origin of magnetic fields in galaxies 0- 
[1] , galaxy clusters [1, 0] , the large-scale structure of the 
Universe [ a an d even the formation of the first stars and 
galaxies |9Hl2|. 

Depending on the environment, the turbulent Mach 
number A^, i.e. the ratio of the turbulent velocity vt 
to the sound speed c^, may vary considerably. In stellar 
interiors and during primordial star formation, we ex- 
pect subsonic turbulence close to the well-studied case 
of Kolmogorov [l3| , while turbulence in galaxies may ex- 
hibit Mach numbers considerably larger than 1, which is 
typical for Burgers-type turbulence jl4| . Observations in 
turbulent molecular clouds revealed the so-called Larson 
laws, hinting towards an intermediate case in between 
Burgers and Kolmogorov turbulence The Larson 

laws consist of: (i) a power relationship between the ve- 
locity dispersion cr^, and the spatial scale of the emit- 
ting volume L, cx L'^ '^*, (ii) selfgravitational cquilib- 
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rium, 2ct„L^/GM cx 1, and (iii) an inverse relationship 
between the mean density n, and size of the cloud, n cx 
L^-^. The Larson laws therefore describe observed cloud 
properties in the turbulent and highly compressible in- 
terstellar medium. 

For the modeling of turbulent magnetic field amplifi- 
cation, the well-studied Kraichnan-Kazantsev model 
[Tsl is often employed. The Kraichnan model [17, 1^ con- 
siders the advection of a passive scalar field, which may 
for instance represent some chemical species, while the 
Kazantsev model [l^ describes the turbulent diffusion of 
a passive vector field. Both models employ the same as- 
sumptions, in particular a velocity field based on a zero 
mean Gaussian random process as well as a ^-correlation 
in time. Both models have been largely used over the 
yea rs and applied to diverse physical circumstances (e.g. 
[19l - l23 | and [l^, [IB] ) and represent an important theoret- 
ical laboratory for studying the multiparticle statistics in 
fluid turbulence. As we are in the following considering 
the amplification of vector quantities, we will adopt the 
term Kazantsev model for brevity. 

The classical studies of the Kazantsev model typically 
focused only on Kolmogorov turbulence. As discussed 
above, the latter is often not applicable in astrophysical 
environments, and steeper slopes for the turbulent spec- 
tra are frequently found both in numerical simulations 
and observational data sets. A systematic exploration 
of such steeper slopes is therefore in order. We pursue 
that here for the full range of magnetic Prandtl numbers 
Pm, denoting the ratio of kinetic viscosity v and mag- 
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netic viscosity 77, from Pm ^ 1 to Pm ^ 1, covering the 
complete range that may occur in astrophysical environ- 
ments. For instance, in stellar interiors, Pm is consider- 
ably smaller than 1, while the gas in galaxies and in the 
intergalactic medium exhibits magnetic Prandtl numbers 
much larger than 1. In primordial clouds, the latter can 
also change during the evolution as a result of ambipo- 
lar diffusion [23|- Over the last years, studies on both 
very large magnetic Prandtl numbers jl, [l9|, [10, [28l - l30| , 
where the amplification is caused by random stretching 
of the magnetic field at the viscous scale, and very small 
magnetic Prandtl numbers [l^, [IJ, l3ll - [33| . where it is 
driven by the inertial-range velocity field, have been re- 
ported. The above-mentioned studies employed a large 
range of analytical and numerical methods, which usu- 
ally are only applicable in a limited regime. Analytical 
studies making use of the Kazantsev model may for in- 
stance treat the regime of Pm ^ 1 and Pm ^ 1, but not 
the regime of Pm ^ 1. Conversely, numerical simulations 
are bound to the regime Pm ~ 1 [341-36] . A comparison 
of these approaches is thus very difficult. As mentioned 
above, most studies so far explored only the effect of Kol- 
mogorov turbulence, although the interstellar medium in 
the first galaxies is expected to be in a highly compress- 
ible regime and potentially closer to Burgers turbulence. 

In this paper, we accurately solve (via numerical inte- 
gration) the Kazantsev equation and calculate the growth 
rate of the small-scale dynamo for the complete range 
of magnetic Prandtl numbers, from the smallest to the 
largest, for different types of turbulence. We define the 
kinematic and magnetic Reynolds numbers as Re=VL / v 
and Rm=VL/?7, where V is the typical velocity at the 
largest scale of the inertial range L. We thus have Pm = 
Rm/Re. We note that, in addition to microphysical pro- 
cesses like Ohmic diffusion, turbulent diffusion process 
may play a substantial role, as pointed out by Weinan 
and Vanden-Eijnden [s?!, Sreenivasan and Schumacher 
[3^ . Here, these effects are incorporated into our model 
in terms of a diffusion coefficient. Covering the range 
sketched above provides solutions for the growth rate 
in the astrophysically relevant regimes, and simultane- 
ously establishes a pathway for a comparison with 3D 
magneto-hydrodynamical simulations. In the following 
Section, we present the Kazantsev model and our nu- 
merical treatment. Subsequently, we show the different 
growing modes for Kolmogorov turbulence, discuss the 
dependence on different types of turbulence, analyze the 
behavior close to the threshold and discuss the implica- 
tions of our results. 



II. KAZANTSEV MODEL AND NUMERICAL 
SOLUTION 

The magnetic field amplification is governed by the 
induction equation 



where B is the magnetic field, and v is the fluid velocity. 

The Kazantsev model PJI] describes the small-scale dy- 
namo process for a Gaussian random velocity field which 
is 5-correlated in time. For this purpose, we introduce 
the covariance tensor as 



((5«, (ri , t)&v, {r2,s)) ^T,,{r)6{t- s) 



(2) 



being Tij{r) the two-point correlation function, defined 
by Batchelor [39| as 



^)T^ir) + ^n{r) (3) 



with Tn^v) and TL{r) being the transversal and longitu- 
dinal part, respectively. We use the model for the correla- 
tion function of the turbulent velocity field has presented 
in Schober et al. 
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l_i?e(i-^)/(i+^) (f)') 0<r<C 

e^<r <L 
L<r, 

(4) 

where £^ = L Re ^/f^'+i) denotes the cutoff scale of the 
turbulence, i. e. the viscous scale, and L the length of the 
largest eddies. 

The transverse correlation function for the general slope 
of the turbulent velocity spectrum is 



' * <r <L 

L < r, 

(5) 

with t{d) = (21 — 38i?)/5. In this formulation we ne- 
glect the effect of the helicity which is not important at 
small scales. We note that, when refering to different 
types of turbulence below, we refer to the velocity cor- 
relation functions defined here, adopting the respective 
slope of the turbulence model, Sv oc These models 
of course only represent an approximation of real astro- 
physical turbulence, where the correlation time is finite. 

By analogy, one can define the components and 
Ml for the correlation function of B, where Ml and M^ 
represent the correlations between the same components 
(e.g. {Bx{0)Bx(r))) and between different components 
(e.g. {Bx{Q)By{r))) of B. In order to ensure the con- 
straint V ■ B = 0, these functions need to obey the fol- 
lowing equation: 



(6) 



SB 
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(1) 



Solutions for the Kazantsev equation have been pre- 
viously proposed in many works and for different pa- 
rameters by assuming isotropy and homogeneity [2^ [2^ . 
In a recent study, Schober et al. [l^ discussed an ana- 
lytical solution based on the quantum-mechanical WKB 
approximation for the case of infinite magnetic Prandtl 
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numbers for different types of turbulence. By using the 
ansatz given in their paper, it is possible to obtain the 
Kazantsev equation from Eq. ([T]) as 

-K(r)— |^ + [/(r)*(r) --r*(r), (7) 

with 'i'{r) related to through the formula 

ML{r,t) = *(r)e2"/(r^V«)- Equation (O formally 
looks like the quantum mechanical Schrodinger equation 
with r being the growth rate and U (r) the potential, 
defined as 

^(,)^^_(^^27^^2(7WV^_ (8) 
^ ' 2 An r ^ ' 

The diffusion of the magnetic correlation K{r) is defined 
as follow: 



Re = 10 



K{r)=rj + TL{0)-TLir). 



(9) 



The latter contains the magnetic diffusivity rj and the 
scale-dependent turbulent diffusion Tl{0) — T^ir) (see 
also Weinan and Vanden-Eijnden [s^l, Sreenivasan and 
Schumacher [sl] for the implications of turbulent diffu- 
sion in the Kraichnan model). It is worth noting that 
numerical simulations show the presence of shear even at 
high Mach numbers, implying that such turbulent diffu- 
sion will be available both in the compressible and the 
incompressible case [13, |4l|. Note also that a solution 
r > for Eq. ([7|) exists only for U{r) sufficiently nega- 
tive in some region of r. The problem is then reduced to 
a search of bound states for U{r). 

Defining x = ln(r) and "^{x) — e^^^9{x), the Kazant- 
sev equation reads 



d^ejx) 

dx^ 



-p{x)e{x) = 



(10) 



with 



p{x) 



(r + U{x))e 
k{x) 



2x 



(11) 



It is useful to introduce the normalised growth rate 



r 



L 



(12) 



with V and L the turbulent velocity and the length of 
the largest eddies, respectively. For further details on 
the derivation of Eq.s (fTO|) and (fTTj) we refer to Scliober 
et al. [3. 

We solve Eq. (|TOl) employing the Numerov algorithm 
[4^ which is well-suited for the study of second order 
problems which contain no first order derivative. Two 
boundary conditions for 9{x) are needed in order to ob- 
tain the solution. We use 

e{x) ^^^^) 0. (13) 

The error in one integration step h is usually 0{h^), 
which leads to a total error in the Numerov method of 
the order 0{h^). We adopt a typical length scale of 1 pc 
(3.18x10^* cm) and a turbulent velocity of 1^ = 1 kms""'^. 
All the parameters have been accurately tested to ensure 
convergence of T and of the wavefunction 9{x). 




FIG. 1: Log-Log plot of the computed fastest and slower 
modes of the normalised growth rate F as a function of Pm 
for incompressible Kolmogorov turbulence and for Re — 10^*. 



III. GROWTH RATE FOR KOLMOGOROV 
TURBULENCE 

In Fig. [T]we report the calculated normalized growth 
rate F for different modes and Re = 10^"* for the case of 
Kolmogorov turbulence from Pm ^ 1 up to Pm 3> 1. 
For a narrow range of Pm we found a strong increase 
of the growth rate, in particular for the fastest growing 
mode, which depends on the fact that, for 5 < Pm < 
10^ the potential is negative both in the inertial and in 
the viscous range yielding two contributions (see Fig. [3]). 
We note that the additional contribution coming from 
the viscous range is marked only for Pm > 10, even if it 
starts to appear for Pm > 5. The higher modes shown 
in the figure clearly depend on the depth of the poten- 
tial U{x), and only the main growing mode (the larger in 
magnitude) was found to exist in the whole range of Pm. 
It is worth noting that the small-scale dynamo amplifi- 
cation can occur also for Pm ^ 1. The presence of the 
higher growing modes is important since it gives an ad- 
ditional contribution to the magnetic field amplification 
that becomes more marked for Pm 3> 1, where the F val- 
ues are 6 orders of magnitude larger than for Pm <C 1 
(for the same Reynolds number). Especially for Pm ^ 1, 
a large number of higher modes has been found. By tak- 
ing in consideration Burgers turbulence we found a small 
number of modes going to a maximum of 3 (Fq, Fi, and 
F2) for Pm —> 00 up to only 1 mode for intermediate 
and small Pm. This provides a further confirmation of 
the fact that for Kolmogorov turbulence we have a larger 
amplification of the magnetic field. 



IV. TURBULENCE DEPENDENCE OF THE 
GROWTH RATE 

In Fig. [5] we report the fastest growing mode for dif- 
ferent types of turbulence indicated by the slope of the 
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FIG. 2: Log-Log plot of the computed fastest normalised 
growth rate F as a function of Pm for five different types 
of turbulence and for Re = 10^''. The viscous and inertial 
regime are explicitely pointed out as well as the intermediate 
range. 

turbulent spectrum Sv{£) oc i'^ , where Sv{£) is the eddy 
velocity at the scale £: Kolmogorov = 1/3), Larson re- 
lation (z9 = 0.38), solenoidal (i9 = 0.43) and compressive 
forcing {■& = 0.48) of turbulence, and compressible Burg- 
ers (i? = 1/2) . Kolmogorov turbulence amplifies the 
magnetic field more effectively compared to the highly 
compressible Burgers case. Furthermore, the growth rate 
linearly increases as Pm^^^ for Pm ^ 1 (for Kolmogorov) 
showing again a flex-point at the boundary of the two 
scales and then becoming constant at very large Pm, con- 
sistent with the WKB results [1^. All the other curves 
show a similar behaviour. 

To support and better explain the results reported in 
Fig. [2]we show in Fig. |3]the negative part of the potential 
U{x) for two different values of Pm and for different mag- 
netic Reynolds numbers. On the upper panel it is clearly 
shown that the potential lies completely in the inertial 
range and the features are led by the magnetic diffusiv- 
ity rj. In the bottom panel the potential in between the 
two regimes is shown. A discontinuity is present at the 
length cutoff £,y, reported as a vertical dashed line, and 
the appearence of a deeper potential in the viscous range 
is the cause of the sudden increase in F for Pm ranging 
from 10^ < Pm < 10^ that we have seen in Fig.s [T] and 
m Again it is important to note that the contribution 
coming from the viscous regime exists also for Pm ^ 5, 
but it is too small and does not affect the integration 
and the final value of F. Only when this contribution be- 
comes deeper than the inertial one, F starts to strongly 
increase. 



V. DYNAMO THRESHOLDS 

An important feature of the small-scale dynamo pro- 
cess is certainly the threshold Rnic, for which the amplifi- 



TABLE I: The normalised growth rate of the small-scale dy- 
namo r as a function of Re and Rm for five different types 
of turbulence. The results for Pm <g 1, and Pm S> 1 are 
reported. 
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cation is activated. We have evaluated Rnic for all types 
of turbulence considered here. An asymptotic value of 
Rmc = 320 has been found for Kolmogorov turbulence 
and Re=10^^, which is in good agreement with other an- 
alytical [3l| and numerical studies [H, |2^, Is^] , which re- 
port a Rmc value of 410, 500, 210, and ~500, respec- 
tively. Furthermore, we found that the threshold for the 
generation of magnetic fluctuations by highly compress- 
ible turbulent flows is considerably larger (Rmc = 32000) 
than for the case of a Kolmogorov fluid. A similar trend 
was also shown in fl^, |4^ . It is worth noting that Leorat 
et al. [44] obtained a critical magnetic Reynolds number 
for compressible Kolmogorov turbulence of the order of 
a few tens. 

In Fig. m we explore the dependence of the fastest 
growing mode on the magnetic Reynolds number for 
values in the vicinity of the threshold Rmc for turbu- 
lence based on Kolmogorov, Burgers and the Larson re- 
lation. We adopt Pm ~ 1 and analyse the scaling of the 
growth rate in the vicinity and far away from the thresh- 
old (Rmc). For all cases considered here, the growth 
rate far from the threshold scales as oc Re" for Pm ^ 1 
(or cx Rm" for Pm ^ 1), with q;=1/2 for Kolmogorov 
and 1/3 for Burgers. In the vicinity of the threshold, the 
growth rate becomes a function of Re- Rmc. Here we per- 
form a fit with a logarithmic function, T — (3 In(Rm) -I- 7, 
with (3 equal to 0.4214, 0.3356, and 0.10438, and 7 = 
/31n(Rmc) to be -2.4516, -2.0901, and -1.0754 for Kol- 
mogorov, Larson and Burgers turbulence, respectively 
[47I . Table U reports the fitted growth rate as a func- 
tion of Re and Rm for magnetic Reynolds numbers far 
from the threshold for the different types of turbulence, 
confirming the most efficient dynamo growth for the case 
of Kolmogorov turbulence. 

VI. CONCLUSIONS 

In this work we computed a numerical solution of the 
Kazantsev equation from the smallest to the largest mag- 
netic Prandtl numbers for different types of turbulence 
with the aim of giving a general view of the full small- 
scale dynamo process. We discussed the features of the 
potential, which determine the behaviour of the growth 
rate. Different types of turbulence, from incompressible 
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FIG. 3; The potential U{x) as a function of the coordinate 
X for two different Pm for Kolmogorov-type turbulence. The 
upper panel shows Pm=10~*, while the bottom panel shows 
Pm=10* and Re=10^^. In both panels the cutoff scale length 
li, is shown as a vertical line. 
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FIG. 4: Computed normalised growth rate F as a function of 
Rm for three different types of turbulence and for Pm ~ 1. 
Fitting curves are also reported in the panels. F scales as 
Rm" (red line) when far from the threshold, and oc /31n(Rm) 
in the vicinity of the threshold. 

Kolmogorov to compressible Burgers, exploring the inter- 
mediate cases of turbulence observed in molecular clouds, 
have been analysed and the threshold for dynamo ac- 
tion been calculated. The growth rate of the magnetic 
energy decreases for more compressible types of turbu- 
lence, consistent with numerical simulations by Federrath 
et al. [4l|- We stress the existence of dynamo action for 
Pm <^ 1 and confirm the earlier results reported in [l^ 
for Pm — > oo. We also included the range Pm ~ 1, where 
the often applied WKB approximation becomes invalid. 
The existence of higher growing modes contributing to 



the amplification of magnetic fields has been analysed 
for the whole range of Pm numbers confirming that a 
larger amplification is expected when it occurs in the 
viscous range. We find that for Kolmogorov turbulence, 
the growth rate scales as Re^/^ for Pm ^ 1, as Rm^/^ 
for Pm ^ 1 and as Ke^^^ and Rm^/'^ for Burgers-type 
turbulence. We further calculated the critical magnetic 
Reynolds numbers Rmc for magnetic field amplification, 
which show a strong increase going from Kolmogorov to 
compressible Burgers turbulence, in good agreement with 
previous analytical and numerical results [24, 29- 31]. 

A caveat that has to be kept in mind is the finite corre- 
lation time that is employed within the Kazantsev model. 
While many of its features, for instance the exponential 
growth and the dependence on the Reynolds numbers are 
in good agreement with numerical simulations, it may 
nevertheless influence its properties to some extend. For 
instance, Schekochihin and Kulsrud [45'] calculated the 
first order corrections of the dynamo growth rate due to 
the finite correlation time, suggesting an overall reduc- 
tion of the growth rate by up to 40%. An additional fea- 
ture that may result from finite correlation times is the 
so-called Golitsyn spectrum [Slj , which appears as an ad- 
ditional term in the correlation function of the magnetic 
field, although it never becomes dominant. Beyond the 
techniques employed in their studies, only the compari- 
son with numerical simulations will shed further light on 
the effect of finite correlation times [33 - 1361 , liol l4l| . 

Our results thus illustrate the behavior of the small- 
scale dynamo under a larger range of astrophysically 
relevant conditions, including the interior of stars, the 
interstellar medium, the intergalactic medium and the 
first stars and galaxies. As originally suggested by [3], 
the initial phase of magnetic field amplification via the 
small-scale dynamo is crucial for providing strong seeds 
on which the a-uj dynamo can subsequently act. Indeed, 
high-resolution numerical simulations find evidence for 
turbulence already in the first protogalaxies, suggesting 
that the dynamo works early on [4^. Due to the ef- 
ficiency of the dynamo even in the highly compressible 
regime, the scenario seems capable of explaining the mag- 
netic field structures in present-day galaxies. 
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